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ABSTRACT 

We use a compilation of cosmic microwave anisotropy data to constrain the epoch of 
reionization in the Universe, as a function of cosmological parameters. We consider 
spatially-flat cosmologies, varying the matter density f2 (the flatness being restored 
by a cosmological constant), the Hubble parameter h and the spectral index n of the 
primordial power spectrum. Our results are quoted both in terms of the maximum 
permitted optical depth to the last-scattering surface, and in terms of the highest 
allowed reionization rcdshift assuming instantaneous reionization. For critical-density 
models, significantly-tilted power spectra are excluded as they cannot fit the current 
data for any amount of reionization, and even scale-invariant models must have an 
optical depth to last scattering of below 0.3. For the currently-favoured low-density 
model with Qq = 0.3 and a cosmological constant, the earliest reionization permitted 
to occur is at around rcdshift 35, which roughly coincides with the highest estimate in 
the literature. We provide general fitting functions for the maximum permitted optical 
depth, as a function of cosmological parameters. We do not consider the inclusion of 
tensor perturbations, but if present they would strengthen the upper limits we quote. 
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1 INTRODUCTION 

The absence of absorption by neutral hydrogen in quasar 
spectra, the Gunn-Peterson effect (Gunn & Peterson 1965; 
see also Steidel & Sargent 1987; Schneider et al. 1991; Webb 
1992; Giallongo et al. 1992,1994), tells us that the Universe 
must have reached a high state of ionization by the red- 
shift of the most distant known quasars, around five. Sev- 
eral mechanisms for reionization, which requires a source 
of ultra-violet photons, have been discussed, and are ex- 
tensively reviewed by Haiman & Knox (1999). In the two 
most popular models, the sources are massive stars in the 
first generation of galaxies, or early generations of quasars, 
and these models have seen quite extensive investigation 
(Couchman & Rees 1986; Shapiro & Giroux 1987; Don- 
ahue & Shull 1991 amongst others). Other possibilities are 
that the reionization is caused by the release of energy from 
a late-decayin g particle, usually thought to be a neutrino 
(Sciama 1993), mechanical heating from supernovae driven 



winds (Schwartz, Ostriker & Yahil 1975; Ikeuchi 1981; Os- 
triker & Cowie 1981) or even by cosmic rays (Ginsburg & 
Ozernoi 1965; Nath & Bierman 199 3) or by radiat ion from 
evaporating primordial black holes (Gibilisco 1996). 

One of the most important consequences of reionization 
is the effect on the anisotropies in the cosmic microwave 



background (CMB), again reviewed by Haiman & Knox 
(1999). Before reionization, the microwave background pho- 
tons have insufficient energy to interact with the atoms, but 
after reionization they can scatter from the liberated elec- 
trons. This leads both to a distortion of the blackbody spec- 
trum and to a damping of the observed anisotropies. Typ- 
ically, the number density of electrons after reionization is 
low enough that only a fraction of the photons are rescat- 
tered, so that some fraction of the original anisotropy is 
preserved. 

There has been continuing rapid progress in observa- 
tions of microwave background anisotropies, and it is now 
well established that there is a rise in the spectrum around 
angular scales of one degree or so, which is where one expects 
to see the first acoustic (or Doppler) peak. While the issue 
of whether or not there is an actual peak, with the spectrum 
falling off again on yet smaller angular scales, remains some- 
what controversial, the existence of significant perturbations 
on the degree scale already indicates that reionization can- 
not have occurred extremely early, as that would have wiped 
out the anisotropy signal. A detailed analysis of the current 
constraints on reionization is our purpose in this paper. The 
earliest analysis of this type was made by de Bernardis et 
al. (1997), and more recently Adams et al. (1998) made a 
specific application to the decaying neutrino model. 
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We will work within the class of generalized cold dark 
matter (CDM) models. We consider a subset, where the dark 
matter is cold and the spatial geometry flat, and we assume 
that the initial perturbations are Gaussian and adiabatic, 
with a power-law form, as predicted by the simplest mod- 
els of inflation. Qualitatively, the COBE DMR detections 
(Smoot et al. 1992; Bennett et al. 1996) provided evidence 
supporting this class of models, by showing that large an- 
gular scale fluctuations have a spectrum close to a scale- 
invariant one. Comparison with a range of observations, in- 
cluding the galaxy cluster number density and the galaxy 
power spectrum, have led to several different recipes aiming 
at concordance, with the CDM models presently providing 
the best framework for understanding the evolution of struc- 
ture in the Universe. 

We allow the possible existence of a cosmological con- 
stant, as supported by recent observations of the magnitude- 
redshift relation for Type la supernovae (Garnavich et 
al. 1998; Perlmutter et al. 1998a,b; Riess et al. 1998; Schmidt 
et al. 1998). We fix the baryon density using nucleosynthe- 
sis (Schramm & Turner 1998). The parameters we vary are 
therefore the matter density Qo, the Hubble parameter h 
and the spectral index n of the density perturbations. In 
this paper, we constrain the amount of reionization as a 
function of these parameters, by carrying out a goodness-of- 
fit test against a compilation of microwave anisotropy data. 
We do not consider the related question of finding the over- 
all best-fitting parameters, and in particular of whether the 
favoured parameter regions are much altered by the inclu- 
sion of reionization, leaving that to future work. 



2 REIONIZATION AND THE OPTICAL 
DEPTH 

2.1 The optical depth 

First we briefly review the relation between reionization red- 
shift and the optical depth. The effect on the microwave 
background anisotropies is mainly determined by the opti- 
cal depth to scattering, and doesn't depend too much on the 
exact reionization history; for illustration we will imagine 
that the Universe makes a rapid transition from neutrality 
to complete ionization. 

If the electron number density is n c , and the Thom- 
son scattering cross-section ot, then the optical depth r is 
defined as 



r(t) 



n c (t) cdt , 



(i) 



where to is the present time. To obtain it as a function of 
redshift, we proceed as follows. First we define the ioniza- 
tion fraction x( z ) = n c /n p , where n p is the proton density. 
Assuming a 24% primordial helium fraction, n p = 0.88 nn, 
where tlb is the baryon number density, the present value of 
which is related to the baryon density parameter by 



8nG 



(2) 



with m p being the proton mass and Ho the Hubble param- 
eter in the usual units. For simplicity, we will assume that 
helium is fully ionized as well as hydrogen; allowing for he- 
lium to be only singly ionized is a small correction (as indeed 




Figure 1. Optical depth for instantaneous reionization at red- 
shift z- lon . From top to bottom the curves are f2o = 0.3, 0.6 and 
1. We took a B h 2 = 0.02 and h = 0.65. 



is allowing for the neutrons at all). Two useful relations are 
the redshift evolution of the electron number density 



n c oc (1 + z) , 

and the time-redshift relation 

They give 
t(z) = n p , (Tt c 



dz' 



{1 + Z ' )2 H{?) X(Z ' ) ' 

where the '0' indicates the present value. 

As long as the dominant matter is non-relativistic 



const ; 



(3) 
(4) 

(5) 
(6) 



h 2 ( z ) 

(1 + zf 
and we can write 

r(z) = r* j* (1 + z' f' 2 \J^X{*') dz' , (7) 
where 

r* = 3Ho " B aTC x 0.88- 0.061 Q B h, (8) 
87rGm p 

the last equality following simply by substituting in for all 
the constants, with the usual definition of the Hubble con- 
stant h. A useful equation for the redshift dependence of 
Q, again assuming only non-relativistic matter and a flat 
spatial geometry, is 

3 



n(z) = n ■ 



(1- 



(9) 



i - n + (i + z) 3 n ' 

For illustration we will assume instantaneous reioniza- 
tion at z — Zi on , so that x( 2 ) = 1 f° r z 5= 2 ion and zero 
otherwise. Equation (TFT) can then be integrated to give 



2r* 

n*ion) = ^r- 



n + n (i 



\3\V2 



(10) 



Sample curves are shown in Figure |l|. Inserting the latest 
permitted reionization redshift, Zi on > 5 from the Gunn- 
Peterson effect, implies only that r exceeds a percent or 
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two for typical cosmological parameters. In order to give an 
optical depth of unity, the epoch of reionization would be at 

-2/3 



100 



/ htt B \ 
Vo.03/ 



n, 



1/3 



(11) 



Therefore, we should expect reionization to occur somewhere 
between 5 < Zi on < 100. 



2.2 Estimates of the reionization epoch 

Estimating the reionization redshift theoretically remains an 
uncertain business. Structure formation in the CDM frame- 
work is hierarchical, with the smallest gravitationally-bound 
systems forming first and the bigger ones appearing later, 
by merging of the smaller structures. When the first fluc- 
tuations enter the non-linear growth regime sometime after 
2 ~ 100 (Peebles 1983), we expect the appearance of the first 
bound objects and therefore, the possible onset of reioniza- 
tion. In most reionization models, the assumed recipe is that 
baryons fall into the potential wells of the developing struc- 
tures in the cold dark matter, forming stars and quasars 
which emit ultraviolet radiation. When this radiation es- 
capes the galaxies, it will ionize and heat the intergalactic 
medium (IGM), and the usual calculations aim to estimate 
when sufficient radiation is available to complete the reion- 
ization. This is already a complex and uncertain calculation, 
made more so if one allows for inhomogeneities which can 
strongly affect the recombination rate (Carr, Bond & Ar- 
nett 1984). Further, we should note that other heating con- 
tributions are not currently excluded (Stebbins & Silk 1986; 
Tegmark, Silk & Blanchard 1994; Tegmark & Silk 1995) and 
may even be necessary. Indeed, it has been claimed from ob- 
servations of the present UV background that it may have 
been insufficient to reionize the IGM (Giroux & Shapiro 
1994), suggesting that collisional heating from supernovae- 
driven winds or cosmic rays could also contribute to early 
reionization. 

Density perturbation growth slows down with time, and 
structures in low-density universes cease growing around 

1 + 2 ~ 1/flo. Therefore, given the present observed mat- 
ter power spectrum, this implies that galaxies formed much 
earlier in low matter density universes. Consequently, reion- 
ization is expected to occur earlier in low-density models 
and, given the bigger look-back time, the optical depth will 
be larger. 

The most extensive theoretical calculations, based on 
the Press-Schechter mass function, tend to show that reion- 
ization occurred after z ~ 50, and that a good guess for 
most models would be 2i on ~ 10 — 40 (Fukugita & Kawasaki 
1994; Tegmark et al. 1994; Liddle & Lyth 1995; Tegmark & 
Silk 1995). Low-density models are towards the top of this 
range and critical-density ones towards the lower end (Lid- 
dle & Lyth 1995). These results have some corroboration 
from numerical simulations (Haiman & Loeb 1997). Specif- 
ically, for ACDM models, Ostriker & Gnedin (1996) and 
Baltz, Gnedin & Silk (1997) show that reionization by pop- 
ulation III stars should have sufficed to reionize the IGM by 

2 ~ 20, although recently Haiman (1998) suggested a lower 
reionization redshift of around 2i on = 9—13 for a flat low- 
density model. If this is indeed the case, then besides the 
determination of 2i on via damping of the CMB anisotropics 
by CMB satellites MAP and Planck, the reionization red- 



shift can be measured directly from the spectra of individual 
sources with the Next Generation Space Telescope (Haiman 
& Loeb 1998) or with 21cm "tomography" with the Giant 
Meterwave Radio Telescope (Madau, Meiksin & Rees 1997). 

In summary, the theoretical uncertainties in estimating 
the reionization redshift are large, and the plausible range 
stretches from just above the Gunn-Peterson limit of z ~ 5 
up to perhaps 40. 



2.3 Spectral distortions from reionization 

In the Thomson limit, where the incident photon energy in 
the electron rest frame is much less than the electron rest 
mass-energy, the blackbody form is preserved by scatter- 
ing. However, the spectrum is measured so accurately that 
one can hope to detect deviations (Zel'dovich & Sunyaev 
1969). The best-known example is the Sunyaev-Zel'dovich 
effect in clusters, which is detectable because of the very 
high electron temperatures in clusters. The reionized inter- 
galactic medium is much cooler, but there is much more of 
it. The amount of distortion of the CMB spectrum is defined 
through the Compton y parameter (Zel'dovich & Sunyaev 
1969; Stebbins & Silk 1986; Bartlett & Stebbins 1991): 



V 



j'-T. - k7; MB ^ 



m c c^ 



n e (7TC dt . 



(12) 



At the epochs of interest, the CMB temperature is negligible 
compared to the electron temperature, which we measure in 
units of 10 Kelvin, denoted T4. If the electron temperature 
is taken as constant, this is the same integral as that giving 
the optical depth, apart from the prefactor. 

With the current limits on this distortion coming from 
the FIRAS experiment, y < 1.5 X 10~ 5 (Fixsen et al. 1996), 
this implies, for typical parameters, 



< 400 T. 



-2/3 ( 



V0.03 



-2/3 



1/3 



(13) 



For the expected typical temperature evolution of the inter- 
galactic medium, it is hard to say much solely from the spec- 
tral distortions about the reionization epoch, except that 
the Universe must have undergone a neutral phase. Steb- 
bins & Silk (1986), Bartlett & Stebbins (1991), Sethi & 
Nath (1997) and more recently Weller, Battye & Albrecht 
(1998) show that almost no reasonable reionized cosmolog- 
ical model violates current spectral distortion constraints. 
Therefore, the information coming from the spatial damping 
of CMB anisotropies, rather than from spectral distortions, 
is crucial to determine the history of the reionization epoch, 
and from now we focus on the anisotropy power spectrum. 



3 THE THEORETICAL MODELS 

As stated in the introduction, our aim is to constrain the 
epoch of reionization for a range of spatially-flat cosmologi- 
cal models. We fix the baryon density at Q b h 2 — 0.02 from 
nucleosynthesis (Schramm & Turner 1998). This is at the 



high end of values currently considered, which makes it a 
conservative choice because decreasing the baryon density 
lowers the acoustic peak and hence permits less reioniza- 
tion. We also do not consider tensor perturbations, which 
contribute predominantly to the low multipoles. As with the 
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baryons, our constraints are conservative in that they would 
strengthen if tensors were included, because including them 
lowers the acoustic peaks relative to the \ow-l plateau. 
The three parameters we vary are 

• Oo in the range (0.2, 1). 

• h in the range (0.5, 0.8). 

• n in the range (0.8, 1.2). 



Our focus is directed towards obtaining upper limits on the 
amount of reionization, though in some parts of parameter 
space there are lower limits too. 

The quantity to be compared with observation is the 
radiation angular power spectrum Ce, which needs to be 
computed for each model. The spherical harmonic index £ 
indicates roughly the angular size probed, 8 ~ l/l. The 
power spectrum is readily calculate d using the CMBFAST pro- 
gram ( Seljak fc Zaldarriaga 1996] ), which allows the input 
of all the parameters we need. However, exploring a multi- 
dimensional parameter space is computationally quite in- 
tensive, and rather than run CMBFAST for every choice of 
the optical depth, it is more efficient and flexible to use an 
analytic approximation to the effect of reionization. This en- 
ables us to quickly and accurately generate Ce spectra for 
arbitrary amounts of reionization. 

The first step is to obtain spectra for the case with no 
reionization, for each combination of our three parameters. 
We take our parameters on a discrete grid of dimensions 
9x7x9. From these spectra without reionization, we can 
generate spectra including reionization using a version of 
t he re ionization damping envelope technique of Hu & White 
(1997). This procedure readily generates accurate enough 



Original last-scattering surface 



spectra for comparison with the current observational data, 
as we will show. However, we do caution the reader that 
this approach will not work once data of improved accuracy 
becomes available. Indeed, as shown by the Herculean 8- 
parameter analysis of Tegmark (1999), the error bars on the 
cosmological parameters coming from the CMB don't seem 
to change very much with the addition of more parameters 
(see also Lineweaver 1998) and in the near future the advent 
of better quality data will make them decrease, introducing 
the necessity for a refined treatment of reionization. 

The underlying physics is the following, illustrated in 
Figure ^. Given an optical depth r, the probability that 
a photon we see originated at the original last-scattering 
surface is exp(— r), the exponential accounting for multi- 
ple scatterings. Those photons will still carry the original 
anisotropy, which we will denote by CJ nt . The remaining 
fraction will have scattered at least once. Their contribution 
to the anisotropy depends on scale. On large angular scales, 
they will have rescattered within the same large region and 
will continue to share the same temperature contrast; this 
is simply a statement that causality prevents large-scale 
anisotropies being removed. On small scales, however, the 
rescattered photons come from many different regions with 
different small-scale temperature contrasts (the small circle 
in Figure ^), and their anisotropy averages to zero. Conse- 
quently we have two limiting behaviours for the observed 
anisotropy Cf s 



cy 



Small 



C? 



exp(-2r)C*l nt Large £. 



(14) 
(15) 



Photons 



Observer 




Figure 2. We see a superposition of photons from the original 
last-scattering surface, and those which scattered. Of the latter, 
those which scattered once originated at decoupling on a smaller 
circle, whose size is given by the time from decoupling to rescat- 
tering. Photons which scatter more than once originate within 
this sphere. 



The factor 2 in the latter is because the power spectrum is 
the square of the temperature anisotropy. 



The reionization damping envelope (Hu & White 1997 



is a fitting function which interpolates between these two 
regimes. For a given r, we obtain the observed spectrum by 



ni c? 



(16) 



where the reionization damping envelope IZe is given in 
t erms of the optica l depth and a characteristic scale £ r by 
flHu & White 1997|) 



1 - exp(-2r) 



(17) 



+ exp(-2r) , 
-0.267, c 2 = 0.581, c 3 = -0.172 



withx = £/(£ r + l) and Ci 
and c 4 = 0.0312. 

The characteristic scale comes from the angular scale 
subtended by the horizon when the photons rescatter 
(i.e. that subtended by the circle in Figure 3). In order to 
obtain a highly accurate result, Hu & White (jl997j) compute 
the characteristic scale £ T via an integral which weights the 
horizon scale with the optical depth, but for our purposes 
the simple formula 



i T = (1 + z ion ) 1/2 (1 + 0.084 In too) 



(18) 



gives sufficient accuracy, where Zi on is given by rearranging 
equation ( |ici| ) . This formula is a fit to the angular size of the 
horizon at reionization (Hu & White 1997, with a sign error 
in their paper corrected). 

We are not quite finished yet, because while the reion- 
ization damping envelope accounts for the loss of anisotropy 
due to scattering, it does not allow for the generation of 
new anisotropies because of the peculiar velocities of the 
rescattering electrons. These create a new, but much less 
prominent, acoustic peak at smaller £ than the original one. 
Because it is a minor feature, it can be modelled simply us- 
ing a Gaussian, the amplitude, width and location of which 
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1000 



Figure 3. Generating a Cp curve including reionization, illus- 
trated for n = 1, h = 0.5, f2o = 1 and an optical depth r = 0.4. 
The top curve shows the spectrum without reionization. Apply- 
ing the reionization damping envelope generates the dotted line, 
and the correction for the new acoustic peak, equation ( |ig| ) , then 
gives the lower solid line. This is to be compared with the exact 
result from CMBFAST for this model, shown as the dashed line. 



depend mildly on the cosmology. The form we choose gives 
the reionized spectrum as 



yobs '-'I 



CODS 
P. — 



In a 



1 - f{t) ' 

where 

m = 4«p(-s5 

A = r(r + 0.16), 
cr = 0.85, 
U = 33f7o + 21ft + 12.5r . 



(19) 

(20) 

(21) 
(22) 
(23) 



The various numerical factors were fits from a comparison 
to CMBFAST output in specific cases. This approach is easily 
accurate enough given the current data, especially as the 
data are given in ST which corresponds to the square root 
of the Ct curve. Figure ^ shows an example compared to an 
exact curve from CMBFAST. 



4 THE OBSERVATIONAL DATA 

In recent years the detection of CMB anisotropics on differ- 
ent angular scales has beco me commonplace. At large scales, 
the COB E measurements (3moot et al. 1992; Bennett et 



al. 1996) constrained the amplitude of the spectrum with 
high accuracy, and to some extent the slope. Since then, 
a plethora of ground-based and balloon-borne experiments 
probing medium and small scales have followed, providing 
increasingly accurate measurements. Although there is still 
quite a large scatter, there is very strong evidence for the 
existence of an acoustic peak at £ of a few hundred, as first 
claimed by Scott & White (1994) and by Hancock & Rocha 
(1997), and therefore a limit on the amount of reionization 
damping which is permitted. 

Our data sample is shown in Figure |i|, and Table 1 lists 
the data points and indicates the sources from which they 
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Figure 4. The observed power spectrum of CMB temperature 
fluctuations. Despite the scatter, there is strong evidence of a rise 
to a peak at an £ of a few hundred. 



were obtained. It is similar to the compilations described by 
Hancock & Rocha (1997) and Lineweaver et al. (1997), and 
several researchers have up-to-date compilations available 
on the World Wide Web. We use updated data and thus our 
sample includes: 

• The 8 uncorrelated COBE DMR points from Tegmark 
& Hamilton (1997). 

• The new calibration of the Saskatoon points (Leitch 
1998); the shared calibration error of these points is small 
enough to be neglected. 

• The new updated QMAP results (de Oliveira-Costa et 
al. 1998). 

• The new OVRO Ring5M result (Leitch et al. 1998). 



data in Table 1 along with the corresponding window 
functions, following the method detailed by Lineweaver et 
al. (1997). In brief, the window functions describe how the 
anisotropies at different I contribute to the observed tem- 
perature anisotropies. For a given theoretical model, they 
enable us to derive a prediction for the ST which that ex- 
periment would see, to be compared with the observations 
in Table 1. 

It has been noted that the use of the \ 2 test can give a 
bias in parameter estimation in favour of permitting a lower 
power spectrum amplitude, as in reality there is a tail to high 
temperature fluctuations. Other methods have been pro- 
posed (Bond, Jaffe & Knox 1998; Bartlett et al. 1999) which 
give good approximations to the true likelihood, though they 
require extra information on each experiment which is not 
yet readily available for the full compilation. We do not use 
these more sophisticated techniques here, but do note that 
as these methods are less forgiving of power spectra with 
too low an amplitude, the results from the analysis give 
conservative constraints on the optical depth. 
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Table 1. The data used in this study, plotted in Figure 
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Figure 5. Contours of the maximum permitted optical depth, as 
a function of n and Ho at fixed Qq. The upper panel shows Co = 1 5 
the lower one Qo = 0-3. Regions to the left of the r max = 
line are excluded, as they do not allow a fit to the observational 
data for any optical depth. The data are more constraining for 
Co = 1, with an upper limit of Ho 65 for a scalar invariant 
power spectrum. 



(I) Tegmark & Hamilton 1997; (2) Ganga et al. 1994; (3) 
Gutierrez et al. 1997; Hancock et al. 1997 (binned); (4) Gundcr- 
sen et al. 1995; (5) Femenia et al. 1998; (6) Tucker et al. 1997; (7) 
Piatt et al. 1997; (8) Piccirillo & Calisse 1993; (9 a ) de Bernardis 
et al. 1994; (9 b ) Masi et al. 1996; (10) Tanaka et al. 1996 (binned); 

(II) dc Olivcira-Costa et al. 1998; (12) Netterfield et al. 1997; (13) 
Wilson et al. 1999; (14°) Scott et al. 1996 and Hancock & Rocha 
1997; (14 6 ) Baker et al. 1998 (15) Leitch et al. 1998. 



5 CONSTRAINTS ON THE REIONIZATION 
EPOCH 

A model is specified by four parameters, fio, h, n and r. 
There is an additional hidden parameter, which is the nor- 
malization of the spectrum. We do not fix this by normaliz- 
ing to COBE alone, but rather seek the normalization which 
gives the best fit to the entire data set. We then examine 
whether each model is a good fit to the data. 

There are A r d a ta = 35 data points. Because we are mea- 
suring absolute goodness-of-fit on a model-by-model basis, 
with one hidden parameter, the appropriate distribution for 
the x 2 statistic has A^ata — 1 degrees of freedom. Nothing 



further is to be subtracted from this to allow for the main 
parameters, as they are not being varied in the fit. To be 
specific, the question we are asking is "If you are interested 
in particular values of Qo> h, n and r for some reason other 
than the CMB data, will the predicted CMB anisotropics 
be an adequate fit to the observations?" . To assess whether 
a model is a good fit to the data, we need the confidence 
levels of this distribution. These are 

X34 < 48.6 95% confidence level ; (24) 

X34 < 56.1 99% confidence level . (25) 

Models which fail these criteria are rejected at the given 
level. We will use the 95 per cent exclusion. Our main focus 
is on limiting reionization, so for each choice of flo, h and n, 
we are interested in the largest value of r, r max , which gives 
an acceptable fit. 

Although we are not concerned with finding the over- 
all best-fitting parameters (which would require variation of 
Qb and ideally the inclusion of tensor perturbations, as in 
Tegmark 1999), we note that the absolute best-fitting model 
in our set, Qq = 0.4, h — 0.6, n — 1.15 and r = 0.3 has a 
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Figure 6. Contours of the maximum permitted optical depth, 
as a function of Oq and Hq with n = 1. 



Figure 7. Limits on the reionization redshift, for f2o = 0.3. 
Reionization must occur late than that indicated by the contour 
levels. This plot assumes instantaneous complete reionization. 



X of 32, agreeing remarkably well with expectations for a 
fit to 35 data points with five adjustable parameters (the 
four mentioned plus the amplitude). These \ 2 values agree 
with other analyses of this type (Lineweaver 1998; Tegmark 
1999), and the best-fit model has parameter values in excel- 
lent agreement with indications from other types of obser- 
vation. 

The upper limits on the optical depth are shown in Fig- 
ures^ and q, for different slices across the parameter space. 
For fio = 1, quite a large amount of otherwise- interesting 
parameter space is now excluded by the CMB data, namely 
the region beyond the r max = contour which will not fit 
the data for any value of the optical depth. For the preferred 
Hubble constant values of around Ho = 65 km s~ 1 , the lower 
limit on n is now around n = 1, severely constraining any 
attempts to salvage critical-density CDM models through 
tilting the primordial spectrum. For critical density with 
n — 1, as commonly employed in mixed dark matter mod- 
els, the optical depth is constrained below 0.3 or even 0.2, 
depending on one's preference for Hq [note that the CMB 
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Figure 8. The same data set of Figure |4|. The plotted curves 
show the best-fit model (fio = 0.4 etc.) and two models that don't 
provide a reasonable fit. For the first model, the high optical depth 
compensates the gain of power at small scales caused by the tilt 
of the spectrum with n > 1. 



anisotropics are hardly altered by introduction of some hot 
dark matter in place of cold (Dodelson, Gates & Stebbint 



1996; )]. 

In the low-density case, the constraints on the optical 
depth are weaker, because the first acoustic peak is pre- 
dicted to be higher in the absence of reionization. However, 
as there is a greater optical depth out to a given redshift 
in low-density models, the constraints on the actual reion- 
ization epoch prove to be quite similar. For fio = 0.3, this 
is shown in Figure |?], which was obtained from the opti- 
cal depth, assuming sufficiently-instantaneous reionization, 
using equation ( |lo| ) . We see that for the most commonly dis- 
cussed n = 1 paradigm, the current limit on the reionization 
redshift is around Zi on = 35, which is just about at the up- 
per limit of the theoretically anticipated range discussed in 
Section 2.2. Future observations may well start to eat into 
that range. 

In Figure ^, we show two typical models which fail to 
fit the data, as well as the absolute best-fitting model. 

As well as describing the results graphically, it is useful 
to having a fitting function for the maximum allowed optical 
depth. A good fit for two particular fio values is given by 
the formulae 

r m ax = 0.03 - (2.9 - 1.5n) (h - 0.65) + 1.9(n - 1) (26) 

[«o = 1] ; 

r m ax = 0.36 - (2.4 - lAn) (h - 0.65) + 1.7(n - 1) (27) 

[fio = 0.3] . 

The second of these is illustrated in Figure []. For general fio , 
a suitable interpolation between these is to interpolate the 
three coefficients linearly in \/fio ( e -g- for the first coefficient 
take 0.76 — 0.73 %/fio and so on). 

There is no simple fitting function for the reionization 
redshift, but an analytic fit is obtained by rearranging equa- 
tion (hd) and putting in the fitting functions for r max . 
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Figure 9. An illustration of the fitting function for Qq = 0.3. 
The lines show, from bottom upwards, n increasing from 0.8 in 
steps of 0.05. The points show the exact results, for the n of the 
line to which they are closest. The worst error on r max is around 
0.02. 



6 SUMMARY 

We have developed an analytic method of generating the Ce 
spectra in reionized models from models without reioniza- 
tion, and confronted models with current observational data 
in order to place upper limits on the optical depth caused 
by reionized electrons. We stress that the constraint is best 
expressed on the optical depth, as the main physical effect 
is that rescattered photons lose their short-scale anisotropy 
and to a good approximation it doesn't matter where the 

grattprincr tool plarp Tn gpnpral t.hp optical Hpptti ra a firnc- 



tion of the complete reionization history, as well as the cos- 
mological model, but at least the first of these dependencies 
can be simplified if it is assumed that reionization happens 
completely and fairly rapidly, in which case the constraint 
can be re-expressed as an upper limit on the reionization 
redshift. 

We considered only a single value of the baryon density, 
at the high end of the preferred range, and did not include 
tensor perturbations. The second of these is definitely con- 
servative, and the first likely to be so, so our results can be 
regarded as rather safe upper limits. However, these quanti- 
ties would in general have to be included if one undertakes 
the more ambitious task of trying to estimate best-fitting pa- 
rameters from the data, rather than delimiting the allowed 
region. Several analyses have been carried out in recent years 
to use available information to constrain the cosmological 
parameters, with the majority neglecting the influence of 
reionization (Ganga, Ratra & Sugiyama 1996; White & Silk 
1996; Bond & Jaffe 1997; Lineweaver et al. 1997; Bartlett et 
al. 1998; Hancock et al. 1997; Lineweaver & Barbosa 1998a, 
1998b). Most closely related to this work are the papers 
of de Bernardis et al. (1997) and more recently Tegmark 
(1998), who investigated how reionization could affect cos- 
mological parameter determination. Our results update and 
extend the former paper, by employing more up-to-date data 
and exploring a wider parameter space. Neither of those pa- 
pers aimed at providing detailed constraints on the epoch of 
reionization, preferring instead to find best-fitting parame- 
ters. Although we have not made a serious attempt at pa- 



rameter estimation, we do concur with those papers that 
the best-fitting models have a blue (n > 1) spectrum and 
significant reionization. 
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